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A simple approximate treatment of statistical equilibrium is developed within the linear a 
sigma model by means of the Hartree factorization technique, providing a simple means 
for sampling initial configurations of the chiral field. These are then subjected to pseudo- 
expansions and their non-equilibrium relaxation towards the normal vacuum is studied. 

ON 

> : 1. INTRODUCTION 

O 

The possibility of forming disoriented chiral condensates in high-energy hadron and heavy-ion 
collisions, such as those anticipated at RHIC, has generated considerable interest in the past 
few years. The underlying assumption is that the collision produces an extended hot region 
| within which approximate chiral symmetry is temporarily restored. The rapid disassembly 

of the system then induces a non-equilibrium relaxation of the chiral field which may lead 
to the formation of large coherent sources of correspondingly soft pions [1-5]. Since these 
disoriented chiral 'domains' have well-defined orientations in isospace the associated pion 
multiplicity distributions display anomalous features. The occurrence of the phenomenon 
may thus provide an observational opportunity for testing our understanding of chiral sym- 
metry. Recent reviews of the topic may be found in refs. [6-8]. 

In order to elucidate the conditions for the occurrence of DCC phenomena and the 
prospects for their experimental detection, it is necessary to carry out dynamical simulations 
of the non-equilibrium evolution experienced by the chiral field as it relaxes from an initially 
very excited state, in which chiral symmetry is approximately restored, towards the normal 
vacuum in which the symmetry is significantly broken. The most popular tool for such 
^ ■ dynamical studies has been the linear a model [9] in which the chiral degrees of freedom are 

described by an 0(4) real classical field with a simple non-linear self- interaction [10-20]. Even 
though this description is relatively simple, ignoring all other degrees of freedom (such as 
those represented by other mesonic species or individual partons), it still presents a significant 
computational challenge. It is therefore of practical interest to develop useful approximate 
methods for solving the equation of motion for the chiral fields and for understanding their 
complicated non-linear dynamics. 

We first describe how it is possible to make a simple approximate treatment of statistical 
equilibrium by employing the Hartree factorization technique. We here largely follow the 
developments made in ref. [21]. Subsequently, following ref. [22], we employ that frame- 
work for understanding the non-equilibrium relaxation of chiral matter subjected to pseudo- 
expansions, with particular emphasis on the conditions for achieving an amplification of the 
low-energy pion modes. A more comprehensive study of DCC observables, and the utility of 
the mean- field approximation for their calculation, is being reported elsewhere [23]. 



2. HARTREE APPROXIMATION IN EQUILIBRIUM 



Most dynamical studies of disoriented chiral condensates have been based on the linear o 
model in which the chiral degrees of freedom are described by the real 0(4) field <f> = (a, tt) 
having the equation of motion 
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0(r,t) = (or(r,t),7r(r,t)) . 



(1) 



The three parameters in the model can be fixed by specifying the pion decay constant, 
f n = 92 MeV, and the meson masses, m n = 138 MeV and m a = 600 MeV, leading to the 
values X= (ml- ml) /Iff = 20.14, v = [{ml - 3ml)/(ml - m 2 )] 1 / 2 ^ = 86.71 MeV, and 
H = mlf n = (120.55 MeV) 3 , with h,c=l [21]. 

As is apparent from eq. flU), the vacuum configuration is aligned with the a direction, 
^vac = (/tj 0)) an d at low temperature the fluctuations represent nearly free a and ix mesons. 
In the other extreme, at temperatures well above v, the field fluctuations are centered near 
zero and approximate 0(4) symmetry prevails. 

In the present discussion, we limit the considerations to macroscopically uniform config- 
urations (chiral matter) and therefore enclose the system in a box with periodic boundary 
conditions. It is then possible to uniquely decompose the chiral field, 



4>(r,t) = 0(t) + 6<j>(r,t) 



4> = < <t>> 



(2) 



The first term, </>, is the spatial average of the chiral field and may be identified with the order 
parameter, while the fluctuations, 5<fi(r), represent elementary quasi-particle excitations 
relative to the constant field. 

By taking the spatial average of the full equation of motion ([[]), it is possible to derive 
an equation of motion for the order parameter [24]. If we subsequently subtract that from 
(|ip and apply a Hartree-type factorization [16,21], we obtain a corresponding equation for 
the field fluctuations. The resulting equations of motion are of mean-field form, 
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= \[ <f>l + ^ $<t> 2 y + 2 ^ 5<j>\ >- 



(3) 



n\ = A[3<^ + ~< S(p 2 y + 2 -< 54>\ y - 



(4) 



□ + n 2 ± 5tf> ± = , n\ = X[ 4>l + ~< 6<p 2 y + 2 ~< 5<j>\ y - v 2 



(5) 



Here 5<f>\\ = 6(f) o <fi is the fluctuation along the order parameter and 5(f) ± is the fluctuation 
perpendicular to 4>. Furthermore, the spatial averages < • > have been replaced by the 
corresponding thermal average -< • >~, evaluated at the given temperature T. We note that 
the effective masses increase with the magnitude of the order parameter 0o as well as with the 
field fluctuations. They are degenerate for 0o=O and vanish at the temperature T c = \[2v. 



Moreover, we always have < < //jj. Since the quasi-particles are thus governed by a 
Klein-Gordon equation, it is simple to obtain the thermal averages self-consistently, 

-<*f>- = sg^/^T- (l = e + "»< (6) 
= 4 = * 2 + ,i. m 

The volume of the box is given by Q and k denotes the wave vector of the individual modes 
in the cavity. The corresponding dispersion relations are indicated as well and the resulting 
effective masses are shown in Fig. 1. 

Eq. (^|) was first derived in ref. [24] and the Hartree treatment is in accordance with 
ref. [16]. Furthermore, we note that the terms < 5(f) 2 > in eqs. (0-§D are sums of the 
field fluctuations in each of the N=4 chiral directions and thus constitute the leading-order 
fluctuation contribution in a 1/N expansion. These are the 'direct' terms that have been 
included in a number of previous DCC treatments or discussions in terms of effective masses 
[10-11,16-19,25-27]. The additional fluctuation terms are the 'exchange' terms and each 
is twice the fluctuation along the particular chiral direction considered (either parallel or 
perpendicular to the order parameter). Although these terms are only of next-to-leading 
order in 1/N, they increase the fluctuation contributions by 2/N = 50% and their effect 
may thus be significant. Recent analyses suggest that the mean-field treatment with all the 
fluctuation terms included is in fact a quite good approximation for both the equilibrium 
properties [21] and for the calculation of DCC observables [23]. 

The statistical properties of chiral matter are most naturally derived on the basis of the 
partition function, 

Z T = J D[<0, 4>}e^ E ^^ = J d 4 ±d A ± W T {±,±) , W T (±,±) ~ q -^ Ko+Vt - TSt) , (8) 

where K = i/}q/2 is the kinetic energy density of the order parameter cf>. The statistical 
weight Wt gives the relative probability for finding the system with a specified value of the 
order parameter <fi and its time derivative i/>. Its simple approximate form contains the 
effective potential energy density Vr, which is shown in Fig. |2|, and the entropy density St 
associated with the quasi-particle degrees of freedom for a given value of 4>q. 

The corresponding free energy density Ft — Vr — TSt governs the distribution of the 
order parameter </>. It depends only on the magnitude 0o an d the disorientation angle 
Xo (the angle between <p and the a direction) and it is easy to calculate in the mean- 
field approximation. Figure ||| shows its appearance along the a axis where its minima are 
situated. For high temperatures, the free energy has its minimum near the origin, a reflection 
of the approximate 0(4) symmetry restoration, and as T is decreased, the minimum moves 
smoothly outwards and finally settles at the vacuum value /„.. 

The location of the minimum in Ft gives the most probable value of the order parameter. 
For the relatively small volumes of relevance in collision experiments, there are significant 



statistical fluctuations around the preferred value and the corresponding thermal distribution 
is determined by the statistical weight Wt in (§p. 
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The effective masses fiii and fi± as func- 
tions of </>o for a range of temperatures, 
calculated in the thermodynamic limit, 
L^oo. For temperatures above T c , the 
curves start at (fro = with degenerate 
values, whereas below T c they only exist 
if is sufficiently large. The correspond- 
ing starting points for /xy are connected 
by the dotted curve and, since flu is then 
nearly independent of T, only the curve 
for T=0 is shown. The locations of the 
vacuum values m a , m n are also shown. 



Figure 1: Effective masses. 




The effective potential energy density 
along the a axis for a range of temper- 
atures. For T<T C the curve starts at a 
minimum value between and v; these 
starting points are connected by the dot- 
ted curve, while the dashed curve shows 
the bare potential Vo for which 54> = 0. 
The potential for other orientations of the 
order parameter can be obtained from 
Vt(0o,Xo)=Vt(0o,O) + H<j) {l - cosxo)- 



Order parameter o (MeV) 



Figure 2: Effective potential. 



The free energy density along the a axis 
for a range of temperatures. The solid 
curves show the results for a number of 
temperatures. For T<T C the curve starts 
at a minimum value of (fio and these start- 
ing points are connected by the dotted 
curve, while the dashed curve shows the 
result obtained when the temperature is 
neglected. The arrow points to the vac- 
uum value i*T=o an d the location of the 
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Figure 3: Free energy. 

The resulting behavior of the magnitude of the order parameter is shown in Fig. f|. As T 
is increased from zero, the fluctuations grow steadily and the equilibrium value of 0o begins 
to decrease from The most rapid change occurs at T « 220 MeV, above which O 
keeps decreasing at an ever slower rate. Correspondingly, increases monotonically with 
T from its free value m n towards ~ 1.6T for T ^> T c , while /iii first decreases, then displays 
a minimum at T « 240 MeV, and finally becomes degenerate with fx± [21]. 




3. SAMPLING OF THERMAL FIELD CONFIGURATIONS 

The simple description resulting from the Hartree factorization makes it possible to de- 
velop a convenient approximate manner for sampling chiral field configurations describing 
macroscopicaliy uniform matter in thermal equilibrium [21]. Since this method is generally 
applicable, it may be of broad interest and we therefore briefly summarize it here. 

The first task is to sample the order parameter [if), <fi) on the basis of the statistical weight 
Wt(iPo, (fio, Xo) given in eq. This quantity factorizes, due to the additive form of the 
exponent. The time derivative ipo is then governed by a four-dimensional normal distribution, 
P^(if)) ~ exp(— VLip1/2T). Furthermore, since the distribution of the magnitude 0o can be 
pretabulated (ignoring at first the H term), the associated sampling task is computationally 
simple. Once <po has been picked, the disorientation angle xo is easy to sample from its 
conditional distribution, P x (xo) ~ exp(— H(j) cos Xo), an d the 0(3) direction ($ , cp Q ) of 7r is 
uniform on Atc. 

Once the magnitude of the order parameter is known, the thermal quasiparticle distri- 
butions are fully determined and the number of quanta in each elementary mode is readily 
sampled, using P(rik) ~ exp(— n-^ek/T) for each of the four principal chiral directions. Since 
the quasi-particle mass tensor is aligned with the 0(4) direction of the order parameter, 
(Xo, $0) fo), a subsequent 0(4) rotation of (f>(r) and if>{r) must then be performed in order 
to express the sampled field configuration in the usual (er, 7r) reference system. 



The most probable value of cfto, the mag- 
nitude of the order parameter, in the 
standard case where H > 0. In the ther- 
modynamics where L —* oo (solid curve) 
4>o is constrained to the value for which 
the free energy density has its minimum 
(see Fig. ^). The bars show the full width 
at half maximum of the thermal distribu- 
tion of 0o in the system with L = 5 fm; 
those for L = 8 fm are about half that 
size. The open dots show the centroids 
for the idealized case having H = 0, for 
the box with L = 5 fm; the behavior is 
qalitatively similar even though the na- 
ture of the phase transition changes. 

Figure 4: Temperature dependence of the order parameter. 
4. EXPANSION DYNAMICS 

We have shown above how thermal equilibrium can be treated approximately in the mean- 
field approximation. However, it is expected that the early collision dynamics causes the 
chiral field to be formed in a state of rapid expansion. The subsequent evolution may 
then lead to a supercooled configuration situated inside the unstable region, thus effectively 
producing a "quench". Several quenched scenarios have been considered [10-12,15-18] but 
they were largely imposed by fiat, thereby reducing the predictive power of the dynamical 
calculations - essentially any degree of magnification can be achieved by suitable adjustment 
of the initial conditions. The degree of arbitrariness can be reduced by examining under 
which conditions a quench-like early scenario may develop dynamically from various possible 
types of initial sceanrios. 

Simple Bjorken-like pictures can be invoked to emulate expansion in D dimensions, ei- 
ther longitudinal (-0=1) [12,17-18], transverse (D=2), or isotropic (-0=3) [13,27]. We have 
studied such scenarios in an approximate manner by augmenting the equation of motion ([[]) 
with a Rayleigh dissipation term [22], 
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if> = d t <t> 



(9) 



The cooling term causes the field fluctuations to decay in the course of time and the associated 
decrease of the energy density is given by E = — (D/t) < tp 2 >. At sufficiently large times, 
the quasiparticle number density decreases as ~ t~ D as t — > oo, as is characteristic of an 
expansion in D dimensions. The time variable should now be reinterpreted as the elapsed 
proper time in a comoving frame (starting at to — 1 fm/c, usually). The corresponding 
Lorentz transformation of the (scaled) spatial coordinates is less essential for our present 
discussion and has therefore been ignored. 



Figures |5]a-5 depict dynamical trajectories for a variety of instructive scenarios. In or- 
der to make a display that does not rely on any assumption with regard to the degree of 
equilibrium, we adopt the field dispersion as a measure of the degree of agitation; it can be 
visualized as a model-independent replacement of the temperature variable. In Fig. is 
shown the dynamical trajectory of the central part of a Ni-sized spherical source prepared 
at T =400 MeV without any initial expansion. The system keeps away from the unstable 
regime, exhibiting an approximately adiabatic evolution. This behavior is rather robust, as 
it occurs for a wide range of initial temperatures and for rod or slab geometries as well. It 
thus appears that initially static field configurations in local equilibrium do not develop any 
instabilities during their subsequent expansion. 

The other trajectories in Fig. [5|a illustrate the effect of endowing the system with an initial 
expansion. The effect increases with D, since the dimensionality of the pseudo-expansion 
effectively acts as the strength of the damping term. The isotropic expansion leads to a 
significant incursion into the unstable region, while the longitudinal one is too slow for that 
and appears to be closer to the self-generated near-adiabatic expansion. 

The approximate equations (|3|-[5D provide a convenient framework for developing an un- 
derstanding of the dynamics generated by the pseudo-expansion (|9]). Imagine that the system 
is initially created in thermal equilibrium at a temperature T well above T c . The field fluc- 
tuations are then sufficiently large to ensure /i 2 >0 in all three eqs. (|]-|). The system is 
expected to experience a cooling resulting from expansion and radiation, so the fluctuations 
decrease in the course of time. This reduces fi 2 which allows the order parameter to grow 
larger, thus counteracting the decrease of the effective masses. The resulting behavior of 
fi 2 is then a delicate balance: for slow cooling the induced growth of 0o is approximately 
adiabatic and the system relaxes towards the vacuum through metastable configurations; 
however, if the fluctuations diminish rapidly a compensating growth of the order parameter 
can no longer occur quickly enough and one or more of the effective masses may turn imag- 
inary, /i 2 <0, indicating that the system has entered an unstable regime where some modes 
grow exponentiallly. 

Figure |5|& shows trajectories for D=3 starting from various temperatures. If the initial 
temperature is lower than 200 MeV or so, the initial value of 0o is already fairly large (over 
60 MeV) and the dynamical trajectories will miss the unstable region. A wide range of higher 
temperatures lead into the unstable region, provided the supercooling occurs sufficiently fast. 
Ultimately, at very high temperatures (above those shown) the system will again stay stable 
throughout, because it takes so long to reduce the fluctuations down to critical size that the 
order parameter has meanwhile had time to start its growth. 




Order parameter o (MeV) Order parameter o (MeV) 

The combined dynamical evolution of the order parameter <pQ ~ < a > and the fluctuations 
< 8(t> 2 > 1 ' 2 . The dashed curve connects the equilibria from T=0 to above 500 MeV and the unstable 
region where fj?^ < is shown by the shaded region; its border intersects the a axis at f n and extends 
up to T c /\/3 at <fto=0. Each system has been prepared in thermal equilibrium at To = 400 MeV, 
using a periodic box (20 fm side length). The irregular solid trajectory (labelled D=0) was obtained 
by solving the standard eq. (|l]) after applying a spherical Saxon- Woods modulation factor (5 fm 
radius and 0.5 fm width) to the hot matter, thereby producing a hot Ni-sizes sphere embedded 
in vacuum; the field was sampled only in the interior (r <2.5). The other three trajectories have 
been obtained by solving the pseudo-expansion equation of motion (|9|) without applying a spatial 
modulation, thus emulating uniform expansions in D = 1,2,3 dimensions. The marks along the 
trajectories are positioned at time intervals of At = 0.2 fm/c. 



Figure 5: Dynamical trajectories. 



T3 



cc 
c 

CO 

CO 
CO 



£ -2 



c 
o 

b. 



i i i i i i i i — ii — | i i i 
\> 7 =400 MeV 




1 2 3 4 5 
Elapsed time t-t (fm/c) 



The time evolution of n\ for pseudo- 
expansions @ with various values of 
the dimensionality D, starting from ther- 
mal equilibrium at To =400 MeV (heavy 
curves), and the corresponding evolution 
of fj,Q (thin curves). The evolution is 
started at to = 1 fm/c, as is commonly 
done, and all curves approach the free 
pion mass (dotted line) for t — > oo. 
Whenever < those transverse modes 
having k 2 < are amplified and the 

resulting maximum degree of amplifica- 
tion is illustrated in Table 1 for a range 
of initial temperatures. 



Figure 6: Time evolution of the pion mass. 

In order to quantify the analysis, it is useful to consider the time evolution of the effective 
masses. Since /ijj > we concentrate on the latter which is illustrated in Fig. |^. It is 
noteworthy that fio ~ A*± throughout the evolution, implying that the amplification of the 
lowest pionic modes is practically identical to that of the order parameter itself. This simple 
feature makes it an easier task to analyze more complicated scenarios. It is convenient to 
express the resulting enhancement of a mode in terms of its amplification coefficient [28], 



G7T 



u>l<0 



dt yJ-U k (t)< 



"1 



k 2 + fi±(ty 



(10) 



which expresses approximately the factor by which the amplitude of a pionic mode has been 
magnified due to incursion(s) into the unstable regime. An upper bound on the magnification 
is provided by the coefficient for k—0, shown in Table [I]. The purely longitudinal expansions 
largely miss the unstable region, while significant magnification occurs for the transverse and 
isotropic expansions, amounting to over a factor of ten in the most favorable cases. 



Table 1: Amplification coefficients/correlation lengths. 



T (MeV) 


D = 1 


D = 2 


D = 3 


200 


0.00 /1.4 


0.02 /1.8 


0.11 /2.0 


220 


0.00 /1.3 


0.50 /1.9 


0.55 /2.5 


240 


0.01 /1.3 


1.20 /2.0 


1.19 /2.7 


300 


0.00 /0.9 


1.84 /1.7 


2.06 /2.7 


400 


0.00 /0.6 


1.67 /1.3 


2.49 /2.1 


500 


0.00 /0.5 


1.31 /l.l 


2.61 /1.6 



The maximum enhancement factor, G£ =0 , in 
macroscopically uniform chiral matter and 
the resulting FWHM of the pion correlation 
function C 7r (ri2) = < 5ir{ri) ■ 5ir{r2) > 
for pions emerging after pseudo-expansions 
@ starting from thermal equilibrium at the 
temperature Tq and using D = 1, 2, 3. 



The relative power spectrum of the pions, 
~ WfcTrjc, where is the Fourier ampli- 
tude of the pion field, plotted as a func- 
tion of the pion kinetic energy uj^-m^. 
The extraction is made at large times 
when the asymptotic scenario of free evo- 
lution has been reached. The plots are 
based on samples of 20 field configura- 
tions prepared at To =400 MeV and sub- 
jected to idealized expansions with either 
D=l (open dots), D=2 (crosses), or D=3 
(solid dots). The irregularities are pri- 
marily due to the shell structure in the 
level density of the cube. 

Figure 7: Pion power spectrum. 

Figure [7| gives an impression of the net effect on the power spectrum of the emerging 
free pions. As expected, the transient instabilities present for D^2 lead to significant en- 
hancements of the power carried off by soft pions. The effect amounts to about an order 
of magnitude for D=3 (relative to the smooth spectrum obtained for D—l), in accordance 
with the amplification coefficients given in Table [T| Although these results were calculated 
for idealized expansion scenarios, they do support the suggestion that such enhancements 
may provide an observable signal of DCC formation [7]. 

5. CONCLUDING REMARKS 

We have discussed some of the key features associated with the formation of disoriented chiral 
condensates in high-energy collisions, such as those planned at RHIC. By application of the 
Hartree factorization technique, it is possible to develop a simple mean-field treatment which 
in turn leads to an efficient method for sampling thermal field configurations. Moreover, the 
mean-field approximation provides a useful conceptual framework for understanding the non- 
equilibrium dynamics of the chiral field as it relaxes from its initial very excited state towards 
the normal vacuum. By augmenting the full equation of motion with a cooling term it is 
possible to emulate chiral matter in uniform Bjorken-like expansion. 

With this method, we have studied the conditions for amplification of the soft pionic 
modes, an important element in the observation of DCCs. Our analysis shows that the oc- 
currence of instabilities, and the associated amplification of pionic modes, depends sensitively 
on the cooling rate, which in turn is intimately related to the character of the expansion. 
Our idealized scenario for D=3 corresponds closely to the isotropic expansion considered 
in refs. [13,27] and our results corroborate the conclusion in [13] that such a scenario leads 
to amplification. Furthermore, our analysis suggest that a longitudinal expansion alone is 
insufficient to cause a quench, if the initial fluctuations are of thermal magnitude. This is 
consistent with what was found in refs. [12,18] for effectively one- dimensional expansions. 
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This qualitative sensitivity to the collision dynamics highlights the importance of employ- 
ing realistic initial conditions for the dynamical simulations of DCC formation. Ultimately, 
the appropriate initial field configurations must be calculated on the basis of the early par- 
tonic evolution, a task which is thus crucial for our ability to assess the prospects of forming 
disoriented chiral condensates in high-energy collisions. 

This work was supported in part by the Director, Office of Energy Research, Office of High 
Energy and Nuclear Physics, Nuclear Physics Division of the U.S. Department of Energy 
under Contract No. DE-AC03-76SF00098. 
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